Giant magnetothermopower of magnon-assisted transport in ferromagnetic tunnel 
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We present a theoretical description of the thermopower due to magnon-assisted tunneling in a 
mesoscopic tunnel junction between two ferromagnetic metals. The thermopower is generated in the 
course of thermal equilibration between two baths of magnons, mediated by electrons. For a junction 
between two ferromagnets with antiparallel polarizations, the ability of magnon-assisted tunneling to 
create thermopower Sap depends on the difference between the size fli,! of the majority and minor- 
ity band Fermi surfaces and it is proportional to a temperature dependent factor [ksT /ud)^^'^ where 
ujd is the magnon Debye energy. The latter factor reflects the fractional change in the net magnetiza- 
tion of the reservoirs due to thermal magnons at temperature T (Bloch's T^/'^ law). In contrast, the 
contribution of magnon-assisted tunneling to the thermopower Sp of a junction with parallel polar- 
izations is negligible. As the relative polarizations of ferromagnetic layers can be manipulated by an 
external magnetic field, a large difference AS = Sap — Sp Sap ^ —{kB/e)f{Ili,'ni){kBT/ujD)^^^ 
results in a magnetothermopower effect. This magnetothermopower effect becomes giant in the 
extreme case of a junction between two half-metallic ferromagnets, AS ~ —kp/e. 
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I. INTRODUCTION 

Spin polarized transport has recently been the sub- 
ject of intense theoretical and experimental interest.^ The 
mismatch of spin currents at the interface between two 
ferromagnetic (F) electrodes with antiparallel spin po- 
larizations produces a larger contact resistance than a 
junction with parallel polarizations, leading to tunneling 
magnetoresistance in F-F junctions^^'^ and giant magne- 
toresistance (GMR) in multilayer structures.'''^ Systems 
displaying GMR have shown other magnetotransport 
effects including substantial magnetothermopower^"^^ 
with a strong temperature dependence. Thermoelectric 
effects have also been discussed in the context of spin in- 
jection across a ferromagnetic-paramagnetic junction.^® 

The Mott formula" S = -(7r2/c|T/3e)((91ncr(e)/(9e),^ 
relates the thermopower 5 of a system to the deriva- 
tive with respect to energy of the electrical conductivity, 
(T(e), near the Fermi energy, ep-, so that, in metals, S 
typically contains a small parameter such as ksT/ep. In 
magnetic multilayers with highly transparent interfaces, 
the Mott formula has been used as a basis for theories 
of transport that explain the origin of the magnetother- 
mopower effect as due to either the difference in the en- 
ergy dependence of the density of states for majority and 
minority spin bands in ferromagnetic layers, ^'^'^^ or a dif- 
ferent efficiency of electron-magnon scattering for carri- 
ers in opposite spin states.* In particular, the electron- 
magnon interaction in a ferromagnetic layer was incorpo- 
rated to explain the observation* of a strong temperature 
dependence of S{T) and gave, theoretically, a much larger 
thermopower in the parallel configuration of multilayers 
with highly transparent interfaces than in the antipar- 



allel one, Sp 3> Sap- For tunnel junctions, magnon- 
assisted processes have been studied both theoretically^^ 
and experimentally^'^ with a view to relate nonlinear I{V) 
characteristics to the density of states of magnons f2((jj) 
as (fl/dV'^ cx fl{eV). 

In this paper we investigate a model of the electron- 
magnon interaction assisted thermopower in a meso- 
scopic size ferromagnet /insulator /ferromagnet tunnel 
junction, which yields a different prediction. The bottle- 
neck of both charge and heat transport lies in a small-area 
tunnel contact between ferromagnetic metals held at dif- 
ferent temperatures. The thermopower is generated in 
the course of thermal equilibration between two baths of 
magnons as mediated by electrons, and, in the relatively 
high resistance antiparallel (AP) configuration of a fer- 
romagnetic tunnel junction, it depends on the difference 
between the size of the majority and minority band Fermi 
surfaces. For a momentum-conserving tunneling model 
we find that 
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e 2n_ 



Sm{T) 



3.47 /keT 
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Sm{T), 



3/2 



(1) 



(2) 



where n_|_ (n_ ) is the area of the maximal cross-section of 
the Fermi surface of majority (minority) electrons in the 
plane parallel to the interface (11+ > n_). The function 
dm{T) is the fractional change in the magnetization of 
the reservoirs due to thermal magnons at temperature T 
(Bloch's r3/2 law), ^ is the spin of localized moments, 
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and lod is the magnon Debye energy. On the other hand, 
we find that the contribution of magnon-assisted tun- 
nehng to the thermopower of a parallel configuration is 
negligible. 

As an extreme example, the magnetothermopower ef- 
fect is most pronounced in the case of half-metallic fer- 
romagnets, where the exchange spin splitting A between 
the majority and minority conduction bands is greater 
than the Fermi energy ep measured from the bottom of 
the majority band, and the Fermi density of states in the 
minority band is zero. In the antiparallel configuration 
of such a junction, where the emission/ absorption of a 
magnon would lift the spin-blockade of electronic trans- 
fer between ferromagnetic metals, we predict a large ther- 
mopower effect, whereas in the lower-resistance parallel 
configuration thermopower is relatively weak:^^ 



the temperature drop, AT. In the linear response regime 
the electrical current may be written as 



Sai 



-0.64^ 

e 



Sp ksT 



Sap 



ep 



(3) 



This is because the contribution of magnon-assisted 
transport to the thermopower in the parallel configura- 
tion Sp is zero and the thermopower only arises from the 
energy dependence of the electronic density of states near 
the Fermi energy. 

The magnon-assisted processes that we consider are 
similar to those discussed in Refs. 19 in relation to the 
magnon contribution to the nonlinear conductance. In a 
ferromagnetic tunnel junction in the antiparallel configu- 
ration, the elastic contribution to the conductance is sup- 
pressed by the mismatch of spin currents at the interface. 
However it is possible to lift spin current blockade while 
conserving the overall spin of the system by emitting or 
absorbing magnons. For example, the change of spin 
that occurs when a minority carrier flips and occupies a 
majority state is compensated by an opposite change of 
spin due to magnon emission. As a result, the spin cur- 
rent carried by electrons crossing an interface between 
oppositely polarized ferromagnets is carried further by 
the flow of magnons (spin waves). 

Microscopically, a typical magnon-assisted process 
that contributes to the thermopower in the antiparal- 
lel configuration, Eqs. (1) and (3), is shown schemati- 
cally in Figure 1. Here the majority electrons on the left 
hand side of the junction are 'spin-up' and the major- 
ity electrons on the right are 'spin-down'. The transi- 
tion begins with a spin- up majority electron on the left, 
that then tunnels through the barrier (without spin flip) 
into an intermediate, virtual state with spin-up minor- 
ity polarization on the right (Figure 1(a)). In the flnal 
step. Figure 1(b), the electron emits a magnon and in- 
corporates itself into a previously unoccupied state in the 
spin-down majority band on the right. In our approach, 
we take into account such inelastic tunneling processes 
that involve magnon emission and absorption on both 
sides of the interface, as well as elastic electron transfer 
processes, in order to obtain a balance equation for the 
current I{V, AT) as a function of bias voltage, V, and of 



/ = GvV + Gat AT, 



(4) 



where Gy is the electrical conductance and Gat is a ther- 
moelectric coeflicient describing the response to a tem- 
perature difference. Under conditions of zero net current, 
the thermopower coeflicient is 



(5) 
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FIG. 1. Schematic of magnon-assisted tunneling via 
an intermediate minority state. This example shows a 
transition from an initial majority state on the left to a 
flnal majority state on the right across a junction in the 
anti-parallel configuration, (a) The process begins with 
a spin-up majority electron on the left, that then tunnels 
through the barrier (without spin flip) into an interme- 
diate, virtual state with spin-up minority polarization on 
the right, (b) The electron emits a magnon (wavy line) 
and incorporates itself into a previously unoccupied state 
in the spin-down majority band on the right. 

The paper is organised as follows. In Section II we 
introduce the model and technique used for describing a 
tunnel junction and we calculate the thermopower in the 
antiparallel (AP) conflguration. We present a detailed 
description of two different models of the interface: a 
uniformly transparent interface where the component of 
momentum parallel to the interface is conserved, and a 
randomly transparent interface. In Section III we demon- 
strate that the contribution of magnon-assisted tunnel- 
ing to the thermopower of a parallel (P) conflguration 
is negligible. In Section IV we discuss the magnetother- 
mopower, give an order of magnitude estimate of the size 
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of the effect, and present results for the magnetother- 
mopower of a junction between two half-metallic ferro- 
magnets. 



II. THERMOPOWER OF AN ANTIPARALLEL 
JUNCTION 

A. Description of the model 

Our initial aim is to write a balance equation for 
the current in terms of occupation numbers of electrons 

«L(R)(eka) = [exp((eka-e^^^V(fci37L(i?.))) + l]"^ and of 
magnons AfL(H)(q) = [exp(Wq/(fcBTi(i{))) - 1]-^ on the 
left (right) hand side of the junction, where 7L(fl) is the 
temperature on the left (right) hand side, e^ — e^ = —eV, 
and ojq is the energy of a magnon of wavevector q. In the 
following we set Tl = T and Tr = T - AT and we shall 
speak throughout in terms of the transfer of electrons 
with charge — e. The index a = {+, — } takes account of 
the splitting of conduction band electrons into majority 
ek+ and minority ek- subbands, eka = efc — aA/2, where 
ek is the bare electron energy and A is the spin splitting 
energy. 

In an AP junction, we assume that the majority elec- 
trons on the left hand side of the junction are 'spin-up' 

and the majority electrons on the right are 'spin-down'. 
The total Hamiltonian of the system is 



H 



h3 



(6) 
(7) 



kk'cK 



where Ht is the tunneling Hamiltonian^^"-^"' that is writ- 
ten in terms of creation and annihilation Fermi operators 
and c. Here a = —a and we assume that spin is con- 
served when an electron tunnels across the interface. The 
tunneling matrix elements ik,k' describe the transfer of 
an electron with wavevector k on the left to the state with 
k' on the right. We will consider ik,k' to be a symmetric 
matrix of the form 



ik,k' = iki 



L2 



1/2 



(8) 



where w|, /{(k) = deL,RO^)/d{hkz) are components of 
electron velocity perpendicular to the interface and 
€L,R(k) denotes the electron energy dispersion in the elec- 
trodes. In our model for t, we neglect its explicit energy 
dependence. However, ik||,k'|| can describe both clean 
and diffusive interfaces by taking into account the con- 
servation of k||, the component of momentum parallel to 
the interface. 

The term Hp^^^ is the Hamiltonian of the ferromag- 
netic electrode on the left (right) side of the junction 
in the absence of tunneling. We use the so called s- 
f (s-d) model,^^'^^ which assumes that magnetism and 



electrical conduction arc caused by different groups of 
electrons that arc coupled via an intra-atomic exchange 
interaction, although we note that the same results, in 
the lowest order of electron-magnon interactions, may be 
obtained from a model of itinerant ferromagnets.^'^ The 
magnetism originates from inner atomic shells (e.g., d or 
f) which have unoccupied electronic orbitals and, there- 
fore, possess magnetic moments whereas the conduction 
is related to electrons with dclocalized wave functions. 
Using the Holstein-Primakoff transformation-^* the oper- 
ators of the localized moments in the interaction Hamil- 
tonian can be expressed via magnon creation and anni- 
hilation operators 6^,6. At low temperatures, where the 
average number of magnons is small < 6^6 ><C 2^ (^ is 
the spin of the localized moments), the Hamiltonian of 



the ferromagnet Hp^^^ can be written as follows 

ttHR) _ rrL{R) , TrL(R) , ttL(R) 



(9) 



Hi^"^^ = ekacLcka, eka = ek - aA/2, (10) 



ka 



The first term He^^\ Eq. (10), deals with conduction 
band electrons which are split into majority ek+ and mi- 
nority ek- subbands due to the s-f (s-d) exchange. The 



-"k-q-l- 



Ck- 



(11) 



(12) 



Hamiltonian H, 



L(R) 



Eq. (11), describes free magnons 



with spectrum ujq which in the general case has a gap 

ujq^o = wq. The third term Hem , Eq. (12), is the 
electron-magnon coupling resulting from the intra-atomic 
exchange interaction between the spins of the conduction 
electrons and the localized moments. 

The calculation is performed using standard second 
order perturbation theory. We write the total Hamil- 
tonian, Eq. (6), a.s H = Hq + V, where the perturbation 
V = Ht + H^^ + H^ is the sum of the tunneling Hamil- 
tonian and the electron-magnon interactions in the elec- 
trodes. First order terms provide an elastic contribution 
to the current that do not involve any change of the spin 
orientation of the itinerant electrons, whilst second order 
terms account for inelastic, magnon-assisted processes. 



B. Elastic contribution to the current 

The first order contribution to the current, in the 

antiparallel configuration, arises from clastic tunneling 
without any spin flip between a majority conduction elec- 
tron state on one side of the junction and a minority state 



3 



on the other. Consider for example an initial state con- 
sisting of an additional majority spin up electron on the 
left with wavevector Icl and energy ekL+- This electron 
can tunnel, with matrix element i^fli ^^^^ ^ minority 
spin up state on the right with wavevector Icr and en- 
ergy eicR-- In addition there is a second process which is 
a transition between the same two states, but in the re- 
verse order, giving a contribution to the current with an 
opposite sign. Together, the two processes give a balance 
equation for the first order contribution to the current be- 
tween the majority band on the left and the minority on 
the right. In addition, there are two first order processes 
that result in transitions between the minority band on 
the left and the majority on the right. Overall, the first 
order contribution to the current is I^p where 

/ + OC 

X S{e -eV - ekRa)|"-L(ekLa) [1 - "fl,(ekRa)] - 
- [1 - 7iL(ekLa)]n-fl(ekRQ)}, (13) 

and a = —a. Neglecting terms that contain the small 
parameter kBT/ep, the current may be written as 



jW «£_y(T+_+T_+). 



(14) 



For convenience we have grouped all the information 
about the nature of the interface into a parameter 7aa' , 

T^c'-Att' E l^k„kR|'^(e-ek,J5(e-ekR„,), (15) 

ktkR 

that is equivalent to the sum over all scattering chan- 
nels, between spin states a on the left and a' on the 
right, of the transmission eigenvalues usually introduced 
in the Landauer formula, ^'^"'^^ although we restrict our- 
selves to the tunneling regime in this paper. Later we 
will employ models of two types of interface explicitly: a 
uniformly transparent interface where the component of 
momentum parallel to the interface is conserved, and a 
randomly transparent interface. 



C. Magnon-assisted contribution to the current 

Below we describe processes which contribute to 
magnon-assisted tunneling. For convenience, we divide 
them into two groups that we label as 'electron' and 
'hole' processes. In 'electron' processes, an increase in 
the number of magnons in one electrode is achieved by 
accepting electrons from the other electrode whereas, in 
'hole' processes, an increase in the number of magnons in 
one electrode is achieved by injecting electrons into the 
other electrode. 







FIG. 2. Schematic of four 'electron' type processes, 
across a junction in the antiparallel configuration, that 
involve transitions from majority initial to majority final 
states via a virtual intermediate state in the minority 
band: (i) and (iii) involve magnon emission on the right 
and left hand sides, respectively, whereas (ii) and (iv) 
involve magnon absorption on the right and left. 

The 'electron' processes are shown schematically in 
Figure 2. The straight lines show the direction of elec- 
tron transfer, whereas the wavy lines denote the emis- 
sion or absorption of magnons. The processes are drawn 
using the rule, appropriate for ferromagnetic electron- 
magnon exchange, that an electron in a minority state 
scatters into a majority state by emitting a magnon. 
The 'electron' processes in Figure 2 involve transitions 
from majority initial states to majority final states via 
an intermediate, virtual state in the minority band. For 
example, process (i), which is the same as the process 
shown in more detail in Figure 1, begins with a spin- 
up majority electron on the left with wavevector Icl and 
energy ekL+- Then, this electron tunnels across the bar- 
rier (without spin flip) to occupy a virtual, intermediate 
state with wavevector kn, in the spin up minority band 
on the right as depicted in the right part of Figure 1 (a) 
with energy CkR-- The energy difference between the 
states is CkR ~ fkL + A so that the matrix element for 
the transition contains an energy in the denominator re- 
lated to the inverse lifetime of the electron in the virtual 
state. For ksT, eV ^ A, when both initial and final elec- 
tron states should be taken close to the Fermi level, only 
long wavelength magnons can be emitted, so that the en- 
ergy deficit in the virtual states can be approximated as 
EkR — CkL + A « A. As noticed in Refs. 25, 33, 34, this 
cancels out the large exchange parameter since the same 
electron-core spin exchange constant appears both in the 
splitting between minority and majority bands and in the 
electron-magnon coupling. 

The second part of the transition is sketched in Fig- 
ure 1(b) where the electron in the virtual minority spin 
up state incorporates itself into a state in the majority 
spin down band on the right by emitting a magnon, which 
is shown as a flip of one of the localized moments. The 
wavevector of the electron in the final state is k' and the 
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total energy of the final (many body) state is ek'+ — Wq. 
Similiar considerations enable us to write down the con- 
tribution to the current from all the 'electron' processes 
in Figure 2. We group the processes into pairs which in- 
volve transitions between the same series of states, but in 
the opposite time order so that they give a current with 
different signs, hence their sum gives a balance equation. 
The contributions to the current of the processes (i) and 



(ii), labelled as / 



(iai) 
AP 



riWn) _ 



-47r2- 
h 



are given by 

f oo I , 



ekL 



xS{e 
x|nL(ekL+) 
- [1 - ?iL(ek. 



kLk'q 

) die ~eV- fk' 



-)]^fl(ek'+)A^fl(q)}, (16) 

where q = Icr — k'. The processes (iii) and (iv) are si- 
miliar to processes (i) and (ii), respectively, except that 
electrons interact with magnons in the left electrode: 



^(iii,iv) 
MP 



+ 00 



de E 

k'kRq 

- Ck 



\tk' 



X 6{e - eV ~ ekR,+ ) 5{e - ek'+ - utq) 
x{ - nR{ek^+) [1 - ?iL(ek'+)] [1 + A^L(q)] + 
+ [1 - ?^i?(ekR+)] nL{ek'+ )iVi(q)}, (17) 

where q = Icl — k'. Energies on the right are shifted by 
eV to take account of the voltage difference across the 
junction. 

An example of a 'hole' process is shown in detail in 
Figure 3 (a) and (b), and Figure 3(c) shows the same pro- 
cess, (v), plus other 'hole' processes (vi), (vii), (viii). The 
'hole' processes involve transitions from minority initial 
to minority final states via a virtual intermediate state in 
the majority band. In contrast to the 'electron' processes, 
an increase in the number of magnons in one electrode is 
achieved by injecting electrons into the other electrode. 
As an example we describe in detail the calculation of 
the matrix element for process (v) shown in Figure 3 (a) 
and (b). The initial state has an additional spin down 
minority electron near the Fermi level on the left (left 
part of Figure 3(a)) with wavevector kL. The first step 
in the transition is the creation of an empty state below 
the Fermi level in the spin down majority band on the 
right with wavevector kp. by the absorption of a magnon, 
wavevector q, to elevate a spin down majority electron up 
to a spin up minority state above the Fermi level on the 
right with wavevector k' (right part of Figure 3(a)). The 
second part of the transition is sketched in Figure 3(b) 
where the spin down minority electron on the left tun- 
nels across the barrier (without spin flip) to occupy the 
empty spin down majority state on the right. The contri- 



butions to the current from processes (v) and (vi), 
and from processes (vii) and (viii), / 



(vii,' 
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AP ' 
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FIG. 3. Schematic of 'hole' type processes of magnon- 
assisted tunneling from an initial minority state to a final 
minority state via an intermediate majority state, (a) A 
typical process begins with the absorption of a magnon 
(wavy line) to elevate a spin-down majority electron be- 
low the Fermi level Ep on the right up to a spin-up mi- 
nority state above Ep, creating an empty state below 
Ep in the spin-down majority band, (b) A spin-down 
minority electron on the left tunnels across the barrier 
(without spin flip) to occupy the empty spin down ma- 
jority state on the right, (c) The same process (v) plus 
remaining 'hole' processes: (v) and (vii) involve magnon 
absorption on the right and left hand sides, respectively, 
whereas (vi) and (viii) involve magnon emission on the 
right and left. 



r(v,vi) _ 
MP — 



+ 00 



dej: 



kL,k'| 



are 



kLk'q 

X 6{e - ekL-) '5(e - eV - ek'- + t^q) 
X |nL(ekL-) [1 - "-p(ek'-)] A^p.(q) - 
- [1 - riL(ekL-)] "p(ek'-) [1 + Nniq)] }. (18) 
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/ + 00 
de 
-co , 



X (5(e- eF- ekR-)'5(e- ek'- + Wq) 
x| - nii(ekR-) [1 - ni(ek'-)]iVz,(q) + 

+ [1 - nfl(eka-)] ni(ek'-) [1 + Nl{cO] }, (19) 

To make our analysis transparent, we rewrite the 
magnon-assisted current as the sum of two parts, 

X _ x(i,ii) , x(iii,iv) , 7-(v,vi) , x(vii,viii) 
^AP—^AP +-^AP +^AP +-^AP 



^AP 

7-spont I rstim 



AP ~ ^AP 

stim 
AP ) 



the first of which, labelled I^p^^, does not contain any 
magnon occupation numbers and represents spontaneous 
emission processes, 



xspont _ e {'^++ + — ) 
~ h 



de / Q{ijj) 
Jq 



X |nL(e) [1 - nR{e - eV - uj)] 

-[l-nLie-uj)]nR{e-eV)}, (20) 

where n{Lu) = '^qS{LU — ujq) is the magnon density of 
states that we assume to be the same on both sides of 
the junction. Since our main aim is to demonstrate the 
existence of an effect, wc choose the simple example of a 
bulk, three-dimensional magnon density of states. We as- 
sume a quadratic magnon dispersion, u)g = Dq^, and ap- 
ply the Debye approximation with a maximum magnon 
energy ljd = D{6tt^ /v)'^'''^ where v is the volume of a 
unit cell. This enables us to express the magnon density 
of states as n{uj) = {3/2)Afuj^/^ /lo^^"^ . The term /J™' 
is only non-zero for finite voltage. Wc calculate it in two 
different limits, small temperature difference AT <C T 
and large temperature difference AT = T,Tit = 0. 

where T{x) is the gamma function and ((a;) is Rieniann's 
zeta function. 

The second term, labelled contains all the 

magnon occupation numbers and it represents absorp- 
tion and stimulated emission processes, 

X { [ni(e) - nnie - eV + w)] [T++Nl{oj) + T__Nr{u;)] 

+ K(e) - riRie - eV - oj)] [T++Nr{u;) + T__Nl{u;)] }. 

(22) 



It vanishes in the limit of zero temperature on both sides 
of the junction, but is non-zero for zero bias voltage in 
the presence of a temperature difference. may be 

written explicitly for arbitrary bias voltage and temper- 
atures. 



rstim ^ ^ 

J 4 P — TTTT ^ 



X {eV{T++ + r__)r(|)C(|) [(fcBTi)3/2+ (fcBTfl)3/2] _ 
- (T++ - T__) r(|)C(|) [ikBnr/'-{kBTnr/']}. (23) 

D. Calculation of the thermopower 

The thermopower Sap is determined by setting the to- 



tal current to zero, 



(1) 

AP 



7-spont 
^AP 



Ifp = 0, and finding 



the voltage V induced by the temperature difference AT, 
Sap = —V/ AT. In general we find 



Sap = 



(24) 



This is the main result of the paper, describing junctions 
between ferromagnets of arbitrary polarization strength 
ranging from weak ferromagnets to half-metals. The fac- 
tors C and B are dependent on the ratio AT/T. We 
evaluate them in the limits of small, AT <C T, and large, 
AT = Tl = T, Tr = 0, temperature difference: 



C 



c(i) 



15/8, 
3/4, 



AT < T 
AT = T 



B = 



3/2, AT<r 
(3-y2)/2, AT = T 



(25) 



(26) 



The function dm{T) in Eq. (24) is the change in the 
magnetization due to thermal magnons at temperature 
T (Bloch's T3/2 law), 36 



1 r°° 

5m{T)^— dujn{uj)NLioj) 
(.J^ Jo 



2C 



r(i)C(i)(^) 



3/2 



(27) 



The thermopower is finite because the current Z^'p" con- 
tains a term (last line in Eq. (23)) that depends on the 
temperature difference. It arises from the difference in 
the thermal distribution of magnons A''^ (w) — NrX^^) and 
the process of thermal equilibration between two baths 
of magnons held at different temperatures, which is me- 
diated by electrons, results in a current. The origin of 
the factor 7^+ — in the numerator of Eq. (24) can 
be understood in the following way. In the 'electron' pro- 
cesses in Figure 2, which contribute to "7++, an increase 
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in the number of magnons in one electrode is achieved 
by accepting electrons from the other electrode. On the 
other hand, in the the 'hole' processes in Figure 3, which 

contribute to 7^ , an increase in the number of magnons 

in one electrode is achieved by injecting electrons into the 

other electrode. Hence the contributions of Tj^+ and T. 

proportional to AT in the current Eq. (23) appear with 
opposite signs. 

The sign of the thermopower, Eq. (24), is specified for 
electron (charge — e) transfer processes and under the 
assumption that the exchange between conduction band 
and core electrons has a ferromagnetic sign. For antifer- 
romagnetic exchange, the overall sign of the thermopower 
would be opposite. For example, processes (i) and (ii) in 
Fig. 2, woTild involve magnons on the opposite side of the 
junction, hence the current I^)^p' would be determined by 
magnon occupation numbers Ni,{c{). Note also that we 
considered a bulk, three-dimensional magnon density of 
states, but in general the magnitude and sign of the ther- 
mopower will depend on the magnon spectrum. 



1. Uniformly transparent interface 

We model different types of interface by introducing 
a dependence of the tunneling matrix elements tkL,kR, 

on the wavevectors Icl = (k^j/j^) and = (k^,fc|j), 



is the component parallel to the interface 



where k^^^j^^ 

and is the perpendicular component. For a uni- 

formly transparent interface of area A such that the par- 
allel component of momentum is conserved, we set the 
dimensionless tunneling factor, Eq. (8), equal to 



= \i\ 1^1,11 1,11 ■ 



(28) 



scattering in the plane parallel to the interface and the 
tunneling matrix element is a matrix element of the to- 
tal scattering potential determined with the use of plane 
waves, 



ki' ,ki^ 



(31) 



where rj is the position of the jth contact with area 
a ~ A|.. A product of two tunneling matrix elements, 
averaged with respect to the position of each defect, will 
be large only if the total phase shift is zero which corre- 
sponds to scattering from the same defect. 



t 



ki;,kHi 



\t\ 



(32) 



where a/ A accounts for a reduced effective area and 



that 



• is an effective transparency. This means 



T 



alia \ ( alia 



(33) 



We assume that the densities of states in the spin band 

a are equal on both sides of the junction so that = 

r_+. In the regime 1 > (n+/n_ -t- n_/n+) Sm{T), the 
thermopower Eq. (24) simplifies as 



-C 



kB (n^ 



TP 



e 2n+n_ 



■ 6m{T). 



(34) 



III. THERMOPOWER OF A PARALLEL 
JUNCTION 



so that 



47r^ \tf ^ min{na, Ua'} 



(29) 



where t represents the transparency of the interface and 
lia is the area of the maximal cross-section of the Fermi 

surface of spins a, 11+ > n_ > 0. Then — 71+ = 

T__ = 47r2 \tf (An_//i2) and T++ = An^ \tf {AU+/h^). 
In the regime 1 (l-|-n+/n_) Sm{T), the thermopower 
Eq. (24) simplifies as 



Sap = —C 



kB (n+ - n_) 

e 2n_ 



5m{T). 



(30) 



2. Diffusive tunnel barrier 

We use a model^^ for describing a strongly nonuniform 

interface which is transparent in a finite number of points 
only, randomly distributed over an area A. Each trans- 
parent point is treated as a defect which causes electron 



For parallel orientation of the magnetic polarizations 
of the ferromagncts, we find that the contribution of 
magnon-assisted tunneling to the thermopower is zero 
(upto the lowest order in the electron-magnon interac- 
tion). We consider the majority electrons on both sides 
of the junction to be spin up and the minority electrons to 
be spin down, and the technical details of the calculation 
of the current are similiar to those described previously 
for the antiparallel orientation. There is a first order, 
elastic contribution / 



(1) 
p ' 



r(i) 



involving tunneling between majority states on the left 
and right, 7^+, and tunneling between minority states, 

T. , without any spin flip processes. The first term in 

Eq. (35) corresponds to a large current response to finite 
voltage whereas the second term arises from the energy 
dependence of the electronic density of states near the 
Fermi energy. 
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(i) 

kL 


k. >^ 


(ii) 

kL 




—t^ 


u ^ 


— — 




(iii) 

k' 


kR 

*4 


(iv) 
k' 


kR 

^4 


^ Lt 


TR 




"R 


(V) 




(vi) 










tR ^ 


(vii) 

k' 


kR 


(viii) 

k' 


kR 







>^ 


^ 

tR 



FIG. 4. Schematic of magnon-assisted tunneling across 
a junction with ferromagnetic electrodes in the parallel 
configuration. Eight processes which, to lowest order in 
the electron-magnon interaction, contribute to magnon- 
assisted tunneling, (i) - (iv) involve transitions between 
minority states on the left and majority states on the 
right via a virtual intermediate state, (v) - (viii) involve 
transitions between majority states on the left and minor- 
ity states on the right via a virtual intermediate state. 

To the lowest order in the electron-magnon interac- 
tion there are eight magnon-assisted tunneling processes 
which are shown schematically in Figure 4. The top 
four processes, (i)-(iv), involve transitions between mi- 
nority states on the left and majority states on the right, 
whereas the lower four processes, (v)-(viii), involve tran- 
sitions between majority states on the left and minority 
states on the right. The overall contribution to the ther- 
mopower is zero because the stimulated emission part of 
the current does not depend on the temperature differ- 
ence across the junction, 

e ?! 



x{ey(r+_ + T_+)r(|)c(|)[(fcsTi)3/2+(fc5r;^,)3/2] _ 

-{r+^-T^^)T{^)C{^)[{ksnf/-'HkBTRf'^]]. (36) 

This can be understood by examining the processes in 
Figure 4. The top four processes (i)-(iv), which have 



minority states on the left and majority on the right, 
all produce terms proportional to T-^. Two of them, 
(i) and (ii), are 'electron' type processes in which an in- 
crease (decrease) in the number of magnons on the right 
is achieved by accepting (injecting) electrons from (into) 
the left, but the other two, (iii) and (iv), are 'hole' type 
processes in which an increase (decrease) in the number 
of magnons on the left is achieved by injecting (accept- 
ing) electrons into (from) the right. Therefore the term 
proportional to 71+ in the last line of Ip™, Eq. (36), 
does not depend on the temperature difference but a sum 
(fcBTi)5/2 _^ {kBTRf/"^. The same is true for the lower 
four processes in Figure 4, (v) - (viii), that produce terms 
proportional to 



IV. CONCLUSION 

As shown above, the thermopower of a tunnel F-F 
junction in the parallel configuration, Sp ~ k^T/ (eep), 
is smaller than the contribution of magnon-assisted trans- 
port to the thermopower Sap in the antiparallel configu- 
ration, Eq. (1). As the relative polarizations of ferromag- 
netic layers can be manipulated by an external magnetic 
field, the large difference AS — Sap — Sp results in a 
magnetothermopower effect. As a rough estimate, we 
take ep = 5eV and Sm = 7.5x 10~^ T'^/^ (for a ferromag- 
net such as Ni, Ref. 36) to give Sap ^ — 3/xVK^^ and 
Sp - 0.5AiVK-i at T = 300K. 

As an extreme example, we predict a giant magne- 
tothermopower for a junction between two half-metallic 
ferromagnets. In a half-metal the splitting A between 
the majority and minority conduction bands is greater 
than ep measured from the bottom of the majority band 
so that only majority carriers are present at the Fermi 
energy. In this case = T_+ = T__ = in Eq. (24) 
and, in the linear regime AT <ti Tl, 



Sap = -0.64- 



Sf 



31, 

eep 



(37) 



This result is independent of temperature and of the spe- 
cific half-metallic material, and it represents a giant mag- 
netothermopower effect AS" « Sap ~ — 55/iVK~^. 

A strong polarization dependence of the thermopower, 
AS* w Sap, enables one to separate the interface con- 
tribution to the thermopower from effects arising from 
a finite temperature gradient in the reservoirs. We as- 
sume in our analysis that phonon-mediated heat conduc- 
tion is much lower than the electronic one, and that a 
fast temperature equilibration inside the ferromagnetic 
metal makes a finite temperature drop across the tun- 
nel barrier possible. Furthermore, the predicted inter- 
face magnetothermopower will be most pronounced in a 
geometry where the bottleneck for electron transport is 
also the bottleneck for thermal transport: in a small area 
mesoscopic junction, ideally, in a suspended STM-type 
geometry. 
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